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' We study the propagation of electromagnetic waves in a spacetime devoid of a metric 

^ " ^ but equipped with a linear electromagnetic spacetime relation H ~ X ' ^- Here H is 

, the electromagnetic excitation {"D./H) and F the field strength {E,B), whereas x (36 

C " 3 ■ independent components) characterizes the electromagnetic permittivity/permeability 

of spacetime. We derive analytically the corresponding Fresnel equation and show that 
it is always quartic in the wave covectors. We study the 'Fresnel tensor density' CJ*J*' as 
(cubic) function of x a-nd identify the leading part of x (20 components) as indispens- 
5h 1 able for light propagation. Upon requiring electric/magnetic reciprocity of the spacetime 

relation, the leading part of x induces the light cone structure of spacetime (9 compo- 
^ ^ nents), i.e., the spacetime metric up to a function. The possible existence of an Abelian 

• i-H , axion field (1 component of x) and/or of a skewon field (15 components) and their effect 

■ on light propagation is discussed in some detail. The newly introduced skewon field is 
' expected to be T-odd and related to dissipation, file nonsym37.tex, 2002-03-24 

1. Introduction 

In pre-metric electrodynamics, tlie axioms of electric charge and of magnetic flux 
conservation manifest themselves in the Maxwell equations for the excitation H = 
{V, H) and the field strength F = {E, B), see I: 

dH ^ J, dF ^0. (1) 

Here J = {p,j) is the electric current. These equations keep their form in all ref- 
erence frames and, since they are metric independent, in all gravitational fields as 
well. In order to develop this scheme into a predictive physical theory, one needs 
additionally a spacetime relation H = H(F), which is a functional in general but is 
assumed here to be of a local nature. It should be well understood, this spacetime 
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H=- H,j dx' Adx\ F = - F,j dx' A dx^ . (2) 
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2 G.F. Rubilar, Yu.N. Obukhov, F.W. Hehl 

relation is a universal constitutive law for the vacuum. We do not treat the consti- 
tutive behavior of material media here, we rather search for a universal law for the 
vacuum, i.e., for the spacetime manifold itself^ 

The simplest assumption is that the spacetime relation is local and linear. If local 
coordinates x^ are given, with ... — 0,1,2,3, we can decompose the excitation 
and field strength 2-forms into their components according to 

-H,,dx'Adx\ F^- 
2 ■' 2 

Then the linear spacetime relation, see i'i, can be written as 

1 

4 

Here iijki is the Levi-Civita symbol, with €0123 = 1- The quantity x^''^^{^) char- 
acterizes the electromagnetic properties of the vacuum and is as such of universal 
importance. It is an untwisted tensor density of weight -1-1 with 36 independent com- 
ponents. If we decompose it into irreducible pieces with respect to the 6-dimensional 
linear group, then we find 

^rjkl ^ (l)^yfci ^_ {2)^ijkl ^ {3)^ijkl ^ ^j^j^ 36 = 20 © 15 © 1 (4) 

independent components, respectively. The irreducible pieces of x ^re defined as 
follows: 

(2}^ijkl ._ 1 f^^ijkl _ ^kli]^^ ^ _(2)^kUj ^ i3)^ijkl ._ ^[ijkl] ^ 

W^ijkl _ ^i]kl _ [2)^ijkl _ {3)^ijkl ^ {l)^klij ^ ^g-j 

Since no metric is available, we cannotform traces. The Abelian axion piece =■ 
a{x) e*^*^' has been introduced by Ni EM. Constitutive laws (for matter) with^^x*"''^' 
(non- vanishing "skewon fields"]^ have been discussed by Nieves and Palllillla. They 
yield P- and CP-violating terms in the field equations. Thus all constitutive func- 
tions in can claim respectability from a physical point of view in the framework 
of a linear response theory. We will come back to this question in Sec. 3. 

Our strategy will now be the following: After having specified the smcetime 
relation in (^), we have a predictive theory. Using Hadamard's method IIj'EJ for 
determining the propagation of waves in the spacetime specified so far, namely a 
4-dimensional differentiable manifold carrying the tensor density x^''''^ of (^), we 
will derive, in Sec. 2, the Fresnel equation controlling the waves. It will turn out 
that the axion a(x) drops out and does not affect the waves, in spite of showing 

''Nevertheless, since Maxwell's equations for material media have the same appearance as the 
vacuum equations (0), a constitutive relation for a material medium can be reminiscent of the 
spacetime relation — and we will profit from this analogy. 

''The namepskewon' has been used earlier in a post-Einsteinian theory of gravity by Mann, Moffat, 
and Taylor see also Soleng and Ee JH. We are using this name in a different context since its 
old meaning doesn't seem to be in use any longer, see a title search in the data banks of SLAC 
and DESY. 
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Generally covariant Fresnel equation and the emergence of the light cone structure 3 

up in the Maxwell equations. In other words, from x^'''^^ in (H)) only W-yAi^^^ the 
leading piece, and (2);;^'^*=' , the skewon field, enter the Fresnel equation. We will give 
(for the first time) a generally covariant and explicit analytical derivation of the 
corresponding Fresnel equation and will show that its wave surfaces are in general 
of quartic order. 

We will study in detail the conditions which must be fulfilled in order to factorize 
the quartic wave surfaces into two equal quadratic wave surfaces, the light cone. 
Up to a function, this amounts to a new derivation of the metric of spacetime from 
electromagnetic data. 

In this article, we will restrict ourselves to linear pre-metric electrodynamics. 
Nevertheless, our methods and results are also useful in media with nonlinear re- 
lations between H and F, provided one studies electromagnetic perturbations on 
top of a background solution fulfilling the field equations. Then one can define, see 
Appendix A, an effective constitutive tensor density x*-^ that is analogous to x*-''^' 
of (|). In this way, our results, in particular our Fresnel equation, also apply to 
the work done in geometrical optics by De Lorenci, Novello, et al.E3 in nonlinear 
electrodynamics a la Born-Infeld, see also 00. 

2. Wave propagation: Fresnel equation 

We will study the pmpagation of a discontinuity of the electromagnetic field follow- 
ing the lines of Ref. Ea. The surface of discontinuity S is defined locally by a function 
$ such that <I> — const, on S. Across S, the geometric Hadamard conditions are 
satisfied: 

[F,j]=0, [d,F,k] = q^ fjk, (6) 
[77y] = 0, [m,k]^q^h,k■ (7) 

Here [T] (x) denotes the discontinuity of a function T across 5", qi := 9^$ is the 
wave covector, and fij and hij are tensors describing the corresponding jumps of the 
derivatives of field strength and excitation. If we use Maxwell's vacuum equations 
dH = and dF = 0, then (||) and (0) yield 

e'^''' qj hki ^ , e'J"'='g,/fe, = 0. (8) 

These equations admit non- vanishing wave covectors q^, provided the constraints 

e''"' f^J fki ^ 0, e'^'''h,,hM=0, e'^''' f,jhki=Q (9) 

are fulfilled. 

We assume that x*"''^' is continuous across S. Then, we find from (^), (|^) and 
®, (0), 

L }_ ^ klmn f (MW 

ij ^ *^ijkl X J mn • V 

With this spacetime relation for the jumps, the conditions (H) can be rewritten as 
X'^'' q, fki = ((^^X^^'"' + ('^X'^'') 9, = , e^'^' q, fki^Q, (11) 



4, 2008 8:32 WSPC/Guidelines nonsym37 



4 G.F. Rubilar, Yu.N. Obukhov, F.W. Hehl 

since (3)-)^-^^' j^^ = q, ^^jki j^^ _ q Ti-^^g axion field a drops out even though 
it enters the Maxwell equations. 

Technically, the system (11) was analyzed in Ref. Ila in the framework of a 



3- (co) vector decomposition of the electromagnetic discontinuity tensor fij. Here 
we present a different, generally covariant derivation of the Fresnel equation by 
applying ideas of Tamm 113. 

As a first step, eq.(pll)2 is solved by 



fi] = qiO-j - qjai . (12) 
The covector is only defined up to a 'gauge' transformation 

fli + Xqi , (13) 
with an arbitrary scalar A. We substitute into (|lT|)i: 

X'"" mkai ^ . (14) 

Because of (|l3|), not all the equations in ( p^ are independent. In order to isolate 
the trivial parts of (p^), it is convenient to pick a specific (anholonomic) coframe 
= a^dx^, here a,P,... = 0,1,2,3: Namely, we identify the zeroth leg of the 
coframe with the wave covector, that is, = g or, in components, e;'^ = qi or 
qa = (1,0,0,0). Then 

X Ifsq^/as — X as — X Oh — . (15) 

'Spatial' anholonomic indices run over a, 6,... — 1,2,3. The zeroth component of 
LSl) vanishes. Thus we find 

flfc = W''''ab = Q, with := x°°°'' . (16) 

The necessary and sufficient condition for the existence of non-trivial solutions 
for ab is the vanishing of the determinant of the 3x3 matrix W: 

W := det W = yabcedefW''''w'"'W^f = ^eafccede/X°"°''x°'°'^X°"°^' ■ (17) 

This yields the Fresnel equation. 

We can rewrite W in a fully 4-dimensional covariant manner. First, we observe 
that the 3-dimcnsional Levi-Civita symbol eabc is related to the 4-dimensional one 
by means of eabc = Coabc- Then we can extend one S-component of the constitutive 
tensors to a fourth summation index. As a result, we find the identity 

OaOd 6&0e OcO/ _ 6/36d 070e 050/ _ 1 Q/36<i 6706 050/ /1o^ 

^OafccA X X — '^aP'ySX X X — 2 ta/375X X X i \^°) 

which holds true due to the (anti)symmctry properties of the Levi-Civita symbol 
and of the constitutive tensor. This allows us to rewrite ( |l7|) as 

^A) — ^ ^ ,«/36(i, ,0706, ,050/ /-iqN 

3! 2 '''' ^'^^^ XXX • (ly) 
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Now we apply the same procedure to the second Levi-Civita symbol and finally 
obtain 

W - ^e^p^s exparX'"'°''x°^^^X°'^^ ■ (20) 
Since e/* = qi, this can, in coordinate components, be written as 

W = J e„ 6,,i„ x""" X'"'' X'^*" q.qmii . (21) 
with 9 det(ei"). We define the fourth order tensor density of weight +1 

G''"' ix) ^ e„ trstu X""'''^ X^l^^" X'^^*" • (22) 

It is totall}j_symmetric, Q^^''^{x) = 0'-'"''^''' (x): ^^id thus has 35 independent compo- 
nents (see c3). With (^), we find the Fresnel equation in the generally covariant 
form 

Q''''\x)moqmi = ()- (23) 



The general Fresnel equation (|2^) is always a quartic equation in qi despite the 
fact that it was derived from a determinant of a_3 x 3 matrix quadratic in the wave 
covectors. This corrects Denisov and Denisov c3 who claim that a particular case 
of the general linear constitutive law may yield a sixth order Fresnel equation; in 
c3, even a Fresnel equation of eighth order is claimed to hold. 

The rest of our paper is devoted to finding conditions that will make the quartic 
Fresnel equation factorize into two quadratic ones and to determine under which 
circumstances the Fresnel equation turns out to be a perfect square. 

For practical calculations, it is convenient to use a 1+3 coordinate decomposition 
similar to that of Refs. Correspondingly, we can rewrite our general result 

dH) as 

W = {QoM + qlqa Af" + qlqaqb M"'' + q^aq^c A/"'" + qaq^cqd A^"''"'') = , (24) 
with 

M a"""" , M" 4^°°°" , A/"'' := 6^?°°"'' , (25) 

j^abc . ^QOabc j^jabcd Qabcd (26) 

In terms of the 3x3 matrices 

C\ \ ^^acd X'""' , T^b" ^ x""^^' e^db , (28) 
the tensor density x*-**^' can be written as a 6 x 6 matrix 

JK_ {Bab Vj 

a 
b 



^ = ( Ca A^" ' ' ^'^^ ^' ^' - ^ ^' ^' ^ • ^^^^ 
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For the Af's, we have expUcitly: 

-hcd {A''" A^^ + A"^ 



_p^(a^6)cp^d _ j[dc^(a^jy^b) ^ ^j[{ab) _^dc „ j^d{a j^b}c'-^ 



M' 



abc de{c\ 



M 



abed 



^ef{c^\gh\d Q 



hf 



(30) 
(31) 



(32) 

(33) 
(34) 



Here the M's, as totally symmetric tensors in 3 dimensions, carry 1©3©6©10© 
15 = 35 independent components. All these results have been verified by using the 
computer algebra system Maple together with the tensor package GrTensoi]^. 



3. Axion and skewon 

Before we draw explicitly conclusions from the Fresnel equation (|2^) or (p^, it 
is instructive to study the structure of Q^i^^ . As we saw, this tensor density has 
35 independent components, the same number as (^);)(-*-''^' + (2)-^^'=' This can be 
understood as follows: Because of (pT|)i, we have 

(35) 



gy'='((3);^)^0. 

This can also be directly verified by substituting ^^^^ 



'X ~ e into (|22|). Accordingly, 
if the spacetime relation (|^) contained only the axion piece ^^^Xj the propagation 
of electromagnetic waves would not be well-behaved. Therefore, in nature, it is 
excluded that x consists merely of ^^^x o;- On the other hand, if we substitute x 
into 5, we find 



g^Okl ^ g^Jkl ((1)^ + (2);^ + 0)^) = g^■Jkl ((1)^ + (2)^) 



(36) 



Note that Q depends on x is a cubic way, that is, our last equation is by no means 
trivial. Consequently, if either ^^^x oi' ^^''x or both pieces exist in nature — and 
this is for sure — then the axion piece doesn't interfere with the local laws of 
the propagation of electromagnetic waves as determined by the Fresnel equation. 
Or expressed in a positive manner: The Abelian axion could be around, the laws 
of electrodynamics are compatible with it and the light propagation wouldn't be 
affected locally. There have been intensive experimental searches for axions, so far 



"See 



http: / /grtensor.org 
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without success However, the axion remains a serious candidate for a particle 
search in experimental high energy physics. 



The axion piece has one further interesting property. The kinetic e 
energy-momentum current in pre-metric electrodynamics reads, see 



ectromagnetic 



''^o.-=\[F A{eo,\H)-H h{ec.\F)] . (37) 

If we substitute the spacetime relation (|^) into it, wc find, in reminiscence of (^6|), 

^S„(x) = ((i)x + (^'x + ^'h) = ((i)x + (^'x) . (38) 

Moreover, '^SqC^'^^x) — 0- Clearly then, the axion docs not contribute to the electro- 
magnetic energy-momentum current. Nevertheless, as Ni 013 has shown, it is possible 
to develop a reasonable theory for the Abelian axion. Accordingly, the axion piece 
^•^^X of the constitutive tensor density x cannot be dismissed a priori. 

Let's turn then to the skewon piece ^^•'x- The name we derived from the skew- 
or antisymmetric 6x6 matrix ^^^x can be mapped to. Commonly, this piece is 
not considered seriously. The conventional argument runs as follows, see Post @. 
Suppose a Lagrangian 4-form L exists for the electromagnetic field. In general 
H ~ dL/dF. If H is assumed to be linear in _F, as is done in (||), then L reads 
L H AF ^ X- F ^ ^^^X -F hF + (^^x -F AF, since the piece with (^^x drops 
out of L because of the antisymmetry (2)x'-'''' = — (2)-,^^'^^ ggg ^ Hence, 

m = L{('^X + ^'h) - H^'^X) + L{^^h) ■ (39) 

The term with ^^-'x eventually becomes the Maxwell Lagrangian, the term with (■^^x 
part of the axion Lagrangian. Since we conventionally assume that all information 
of a physical system is coded into its Lagrangian, we reject (2)-)^«ife' ^ as being 
unphysical. 

However, as we saw already in (^8|), if a piece drops out from a certain expression, 
it does not necessarily imply that this piece lost its right of existence. Eq.(^9|) only 
shows that L is 'insensitive' to the skewon piece. In contrast, the electromagnetic 
energy-momentum current ( |3^ ) definitely 'feels' the contribution from the skewon, 

'^I]„((2)x)y^0, (40) 

that is, the skewon piece does carry electromagnetic energy-momentum. Thereby, 
it could displays its presence. Clearly then, we expect that ^^^x influences light 
propagation. And, indeed, it does. By inspection, we find 

g^^^i (d)^ + (2)^) = gr^M ((1)^) ^ ^^^i ((1)^^ (2)^) (41) 

with 

r''(«x/^^x)v^O, (42) 



even though 



g''^^i}'^\) = (). (43) 
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In other words, like the axion, the skewon cannot propagate hght decently unless 
the '-^^x piece participates. Moreover, 

^*jfei((2)^^(3)^) (44) 

i.e., the skewon and the axion piece cannot exist alone or together, the 'lead- 
ing' f^^^x piece is indispensable. If we knew more about the detailed structure of 
g^-' {^^^Xi ^^■'x)j then we could better judge which modes of the skewon are com- 
patible with present-day experiments. 

In any case, in linear pre-metric electrodynamics a Lagrangian is not assumed 
to exist a priori. Therefore, it is not alarming that ^^^x drops out from the proto- 
Lagrangian L. Since it carries electromagnetic energy-momentum, we will take the 
possible existence o/^^^x*-''^' for granted. 

What is then the possihlfi jihysical meaning of '^■'x? This has already been 
discussed by Nieves and Pal lll'tHl. But let us first collect some formalism. Consider 
a certain vector field ^ = ^"Bq, with the basis Cq of the tangent vector space at 
each point of spacetime. Then we can transvect the energy-momentum current ( ^7| ) 
with 

Q := r = i [F A (CJi?) - A (^JF)] . (45) 

The scalar- valued 3-form Q is expected to be related to conserved quantities pro- 
vided we can find suitable (Killing type) vector fields ^. Therefore we determine its 
exterior derivative and find after some algebra, 

dQ^{^\F)hJ+]^{F hC^H-H AC^F) , (46) 

or, in holonomic components, with e^-''^''Qjki/G, :— e^^^^ Jjki/&: and T-C^ := 

a. = ^"Fki J' + ^ {Fki C^n"' - H"' C^Fki) . (47) 

Here £5 denotes the Lie derivative along the vector ^. Now we substitute the linear 
spacetime relation (^), or Ti^^ = x''''™"F„j„/2, and find 

d^Q' = ^Fki + \ [Fki /:5(x'*""^;»0 - x''™"^m„ C^Fm] . (48) 
o 

We apply the Leibniz rule of the Lie derivative and rearrange a bit: 

^^Q' = ^Fm J' + \ [{Ca"'''') FjFm + ix^'"' - x''''')Fj C^Fki] . (49) 
o 

We substitute the irreducible pieces of x*-''^'. Then we have 

= eFki J' + l^e {^'h''"' + ^'^X'^'"') F,,Fki + i ^^h''"' F,, C^Fki . (50) 

Here it is manifest that '■^•'x and the axion ^'^^x behave qualitatively different as 
compared to the skewon ^^■'x- If ''^^X and ^^•'x carry a symmetry along ^ such that 
£^(^^X*"''^' = and £^(3)T^iifci _ then, in vacuum, i.e., for — 0, we still have 
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Generally covariant Fresnel equation and the emergence of the light cone structure 9 

non-conservation because of the offending term ^^^xFF. Here the dot symbohzes 
the Lie derivative along ^. If ^ can be interpreted as a 'time' direction, then Q 
represents the electromagnetic energy density. 

In any case we see that ^^•'x induces a dissipative term with a first 'time' deriva- 
tive. This is what we might have expected since, in general, dissipative phenomena 
cannot be described conveniently in a Lagrangian framework. It is then our hy- 
pothesis that the skewon piece ^^^x can represent a field which is odd imdEr T 
transformations. This is also what had been discussed by Nieves and Pal O'tHl. Of 
course, we must investigate how this skewon^ as we may call it in a preliminary 
way, disturbs the light cone and whether there is perhaps only a viable subclass 
of the 15 independent components of the skewon. An isotropic skewon s could_he 
constructed by putting in (^) C;, = s^t ^a'' — —s5\, see Nieves and Palll3. 

Let us try then to collect our results in order to get a rough picture of the phys- 
ical meaning of these different irreducible pieces. The piece (i)^'-'''' is indispensable 
for an appropriate propagation of light, as we saw in (|35|), (^3|), and (Q). Thus we 
call it the leading piece of x*-''^'. 

If (i)x*^'^' exists alone and for the spacetime relation (^) electric/magnetic reci- 
procity is assumed to hold additionally, see Sec. 4, then, up to an unknown f unction 
the metric of spacetime can be derived including its Lorentzian signature EjqEj. 
One can think of this reduction in the way that electric/magnetic reciprocity cuts 
the 20 components of (i)x*"''^' into half, that is, only 10 components are left for 
the metric. Modulo an undetermined function, we have then 9 remaining compo- 
nents. These 9 remaining components of (i)x*'''^' determine the light cone at each 
point of spacetime. Accordingly, in (i)x*"''^' the light cone of spacetime is hidden 
and thereby conventional Maxwell-Lorentzian vacuum electrodynamics as well. To 
put it more geometrically, the first irreducible piece '^-'xi "^i^- electric/magnetic reci- 
procity, yields the conformal structure of spacetime. In this sense, there is no doubt 
that (I'x*-''^' is the principal part of the constitutive tensor density x of the vacuum. 
But, as we argued above, there are no a priori reasons for excluding either the axion 
(3)-^ or the skewon ^^■'x- In particular, the Fresnel equation and the electromagnetic 
energy-momentum current can accommodate both pieces in a plausible way. 

4. Almost complex structure 

Let us now search for a condition restricting the constitutive tensor density j^i^^ of 
the spacetime relation (H). Since the times of Maxwell and Heaviside, in the equa- 
tions of electrodynamics a certain symmetry was noticed between the electric and 
the magnetic quantities and was used in theoretical discussions. We formulate elec- 
tric/magnetic reciprocity as follows i: Given the energy- momentum current (^7|). 
It is electric-magnetic reciprocal, i.e., it remains invariant ''Sq, — > under the 
transformation 

H^QF, F^--H, (51) 
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with the dimensionful pseudo-scalar function Q — C,{x). 

We now postulate electric/magnetic reciprocity of the spacetime relation 
Under (|5l|), the spacetime relation transforms into 



1 

4C 

If we substitute this into (||), we find, after some algebra, the consistency condition 
1 

with the function 

1 

96 

In order to allow for a more compact notation, we introduce the dimensionless 
density x*^''' := x^'^'^VC and define 

kl p ^mnkl 



C Fij — ^TT ^ijkl x'^'™" Hjnn ■ (52) 



^ iir,rr.n X""^^) (ip.rs X'^'') = , (53) 



C • — n/:? {^ijmn X i^pqrs X ■ C^'^) 



2 

Then the closure relation reads 



X™""' . (55) 



l^ij l^mn — '^^ij (56) 

or, even more compactly, as 6 x 6 matrix equation, 

KK.= 1q. (57) 

Mathematically this means that the operator k represents an almost complex struc- 
ture on the space of 2-forms. For the snccial case of ^'^^ = '^^•'x = 0, such closure 
relations were first discussed by Toupin E2I, Schonberg E3, and Jadczyk Ea. 

Let us now make the closure relation explicit. We turn back to the constitutive 

l— — I o 

6x6 matrix (|29|). We define dimensionless 3x3 matrices A '■= A/C, etc. In terms of 
these dimensionless matrices (we immediately drop the small circle for convenience), 
the closure relation reads, 

A'"'B,b + C\C%^-Sl (58) 
CcA'''' + A^'Vj" ^ 0, (59) 

BacC^b+'Da^B^t^O, (60) 

BacA'^' + V.^Vj' = ~dl (61) 

Assume det B 0. Then we can find the general non-degenerate solution. We define 
the matrix Kat by 

K:^BC, i.e. C = B'^K , (62) 
and substitute it into (|60|), 

V = -KB~\ (63) 
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Generally covariant Fresnel equation and the emergence of the light cone structure 11 

Next, we solve (^8|) with respect to A: 

A^-B-^ -B-^KB-^KB-^. (64) 

We multiply ( |64[) b y C from the left and by T) from the right, respectively, and find 
with (H) and^^, 

CA = -B-^KB-^ - B-^KB-^KB-^KB-^, (65) 
AV = +B-^KB-^ + B-^KB-^KB-^KB-^. (66) 



Thus, (|59| ) is automatically satisfied. Accordingly, only (61) has still to be checked. 
We compute its first and second term of its left side, 

BA^ -I- KB-^KB-^ , (67) 
V"^ = KB-^KB-^ , (68) 

and find that it is fulfilled, indeed. 

Summing up, we have derived the general solution of the closure relation ( [57| ) 
in terms of two arbitrary matrices B and K as 

A = -B-^ - B-^KB-^KB-^ , (69) 
C = B-^K, (70) 
V = -KB-'^. (71) 

The solution thus has 2 x 9 = 18 independent components. Alternatively, one can 
write the solution as 

A = -{l + C'^)B-\ (72) 
V = -BCB-^ , (73) 

which is parametrized by the arbitrary matrices B and C with altogether 18 inde- 
pendent components. 



5. Spacetime metric 

Can we construct the spacetime metric by using the new general solution for the 
generalized closure relation? Technically, this is equivalent to the question: Can the 
general Fresnel equation be reduced to the light-cone structure? Here we study the 
conditions for such a reduction. 

Let us decompose the arbitrary matrix B into its symmetric and antisymmetric 
parts, 

Bab = bab + ^abcU" , with bah ■= B(ab) , := e""-^ B[ab] ■ (74) 

Note that bab contributes to '^^^x and n'^ to ^^^x- Now we can lower the index of 
n° by means of bab, namely, Ua := bbc'n'^ and -n? := rfric = babn'^n'', and, provided 
det6 7^ 0, we can raise an index by 6"'' which denotes the inverse of bab- We find 
det B = det b + n'^, and the inverse of (|7j) reads 

B"'' = — ^ ^ (5"'' + n^n'' - e°^^ n,) . (75) 

det b + ^ ' 
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Here the symmetric matrix 6°^ is the matrix of the minors of hab- If det 6 7^ 0, then 

lab ^ ^ab ^ 

Let us find out what is quahtatively new in the asymmetric case as compared 
to the previously studied symmetric case For this purpose we consider the par- 
ticular solution of the closure relation for = 0. Consequently, C = V = and 
A = -B^^. Then, by substituting (|7|) into @), we obtain: 



ql - 2ql (g2 det h - {qnf) + {q^ det b + {qnfy 



(76) 



(det6-|-n2) 

Here we used the abbreviation (qn) := qan"-- 

In general, this expression is neither a square of a quadratic polynomial nor a 
product of two quadratic polynomials. In other words, neither a light cone nor a 
birefringence (double light cone) structure arises generically. In order to study the 
reduction conditions, let us assume that the Fresnel equation is a product of two 
quadratic equations for qi, i.e., the spacetime 'medium' is birefringent. Accordingly, 
for ([76|) we make the general ansatz 



W 



Ml 



«)('7o-/3) 



[qt-{a + (3)ql+af3] , (77) 



(det& + n2)'^" (det& + n2) 

with some polynomials a and (3 of order 2 in g^. This implies the relations 

a + (3^2{q^~ {qnf) , a/3 = {cf + {qnff , 

with (f := qaQbb"'''- Since a and /3 enter symmetrically in ([ ' 
nonlinear system can be given in the form 

2 



{qn) 



(3 



(qn) 



(78) 

the solutions of this 
(79) 



Thus, the question of the reducibility of the Fresnel equation translates into the 
algebraic problem of whether the square root \/ —if is a real linear polynomial in 
qa- There are three cases, depending on the rank of the 3x3 matrix bah- 

(i) When bab has rank 3, in other words, when det 6 ^ 0, then we can write 
9^ — QaQb b"''' det b, and the general conclusion is that no factorization into light- 
cones is possible (the roots a are complex), unless n° = 0. This latteii_nQndition 
implies that the constitutive tensor is symmetric, and the results of tJ'BEa are 
recovered. 

(ii) When bab has rank 2, i.e., det 6 — 0, but at least one of the minors is 
nontrivial. Then, without loss of generality, we can assume the following structure 
of the matrix b: 

'bn 612 0\ 

bab = \ 612 622 . (80) 
0/ 



Its only non- vanishing minor is 

r33 I 611 612 

b21 &22 



^11^2 



&?2 ^ 0. 



(81) 
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Note that (|8^) is the most general form of a rank 2 matrix hab, up to a renaming of 
the coordinates. In order to avoid complex solutions, we have to assume that the 
minor b — — < 0, so that \/—cp = iiq^. Then ( [79|) leads to 

a = - [qm^ + q2n^ + qsin^ + A^)] ^ = " [^i"^ + 92?^^ + qain^ " A^)] ^ • (82) 

The interpretation is clear: we have birefringence, i.e., two light-cones. In this case, 
the Fresnel equation is found to be 

2 

Oll(?^) +20i2n^n'^ + 022("-) 

X (go + + 92"' + qzin"^ - ') 4^3) 

Then we can read off, up to conformal factors, the components of the two corre- 
sponding 'metric' tensors defining the light-cones: 













\ 









1 9 

n n 







1 9 




(n2)2 




Vo 


'n}{in? + 






{n^ + / 


/I 










\ 









1 9 

n n 


n^{n^ — fi) 





1 2 

n n 




(n2)2 


n n 


Vo 


n^{n^ — 









(84) 



(85) 



We can verify that det(gi-') = det(5^^) = (ni)2(n2)2&33 ^ _^„i)2^^2)2^2 ^ g, so 
that both metrics have the correct Lorentzian signature. 

(iii) When the 3x3 matrix bat has rank 1: In this case all the minors are zero, 
i.e., q""^ — 0, which corresponds to the case 2 for /i 0. We then see that the Fresnel 
equation reduces to a single light cone, but the resulting metric is degenerated, since 
det(5*^) = 0. 



6. Discussion and conclusion 

In this paper we extended appreciably our earlier results of linear pre- 

metric electrodynamics by relaxing the symmetry and the closure property of the 
constitutive tensor x^-''^^ of spacetime. In the most general case of a linear spacetime 
relation (when both closure and symmetry are absent), the wave propagation is 
governed by the newly derived Fresnel type equation ( |2^ ) which is still of quartic 
order, contrary to the claims in Rcf. EJ. We have deduced the tensor density Q'^i'^^ 
in ( p2|) which induces the generally covariant form of the Fresnel equation ( ^3|) . We 
studied some properties of C^'^'^' and showed that (i);^^'^', as the leading part of 
-^tjki ^ is indispensable for a decent propagation of electromagnetic waves. However, 
also each of the two remaining parts, the axion (3)-^yfe' g^jj^j^ ^j^g skewon (2)T^«jfei^ 
has a well-defined physical meaning. The axion drops out from light propagation. 
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High energy physicist search for such type of particles. The skewon seems to be 
T-odd and related to dissipative processes. Its influence on light propagation, see 
(HI), deserves further study. 

In Sec.^, the general solution of the closure relation is presented for the case of 
an asymmetric linear Within a particular class of the asymmetric solutions 

of the closure relation, the reduction of the Fresnel equation was analyzed in detail. 
We have demonstrated that some of these solutions can yield birefringence. A pre- 
liminary study of the general case with K ^ Q shows that the conclusions remain 
qualitatively the same as those presented in Sec]|. 

Thus, we can conclude that the conditions of closure and symmetry of -^^i'^'' 
are sufficient for the existence of a well-defined light cone structure. If any of these 
conditions is violated, the light cone structure seems to be lost. The necessary 
conditions have still to be found. Of course, if one wants to study possible violations 
of Lorentz invariance by means of the skewon (2);^'^'=', e.g., then the light cone 
structure cannot be considered as sacrosanct any longer. 
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Appendix A. An effective constitutive tensor 



Consider an arbitrary local spacetime relation H 
corresponding Maxwell equation, for J — Q, reads 



kl 



0. 



H{F). In components, the 



(A.l) 



A small perturbation AF of the electromagnetic field around some background F 
can be written as -F = -I- AF. Then, to first order in the perturbation, we have 
for the excitation 

1 dHki 



2 dF„ 



AF„ 



(A.2) 



Inserting ( [A. 2D into (A.l) and assuming that the background field F is a solution 
of (A.l), i.e. e'-''^' dj Hki{F) = 0, we obtain an equation for the perturbation: 



ijkl 
Xcff 



ijnin 



dH„ 



dFu 



(A.3) 



The effective constitutive tensor density x^s^ will, in general, depend on the local 
constitutive law and on the background field F . 
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